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Abstract 

The 1 = 2 -kit scattering length is calculated in fully-dynamical lattice QCD with domain-wall 
C > . valence quarks on the asqtad-improved coarse MILC configurations (with fourth-rooted staggered 

sea quarks) at four light-quark masses. Two- and three-flavor mixed-action chiral perturbation 
theory at next-to-leading order is used to perform the chiral and continuum extrapolations. At the 
^ physical charged pion mass, we find m^a 1 ^ 2 = —0.04330 ± 0.00042, where the error bar combines 

■ the statistical and systematic uncertainties in quadrature. 
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I. INTRODUCTION 



Pion-pion (tttt) scattering at low energies is the simplest and best-understood hadron-hadron 
scattering process. Its simplicity and tractability follow from the fact that the pions are 
identified as the pseudo-Goldstone bosons associated with the spontaneous breaking of the 
approximate chiral symmetry of Quantum Chromodynamics (QCD). For this reason, the 
low-momentum interactions of pions are strongly constrained by the approximate chiral 
symmetries, more so than other hadrons. The scattering lengths for tttt scattering in the 
s-wave are uniquely predicted at leading order (LO) in chiral perturbation theory (%-PT) 

my=° = 0.1588 ; m^a{f = -0.04537 , (1) 

at the charged pion mass. Subleading orders in the chiral expansion of the tttt amplitude 
give rise to perturbatively-small deviations from the tree level, and contain both calculable 
non-analytic contributions and analytic terms with new coefficients that are not determined 
by chiral symmetry alone @, 0, H[. In order to have predictive power at subleading orders, 
these coefficients must be obtained from experiment or computed with lattice QCD. 

Recent experimental efforts have been made to compute the s-wave tttt scattering lengths, 
a%° (/ = 0) and a{= 2 (1 = 2): E865 @, Q (K e4 decays), CERN DIRAC Q (pionium 
lifetime) and CERN NA48/2 {K* — > 7T ± 7r 7r°). Unfortunately, these experiments do not 
provide stringent constraints on a 1 ^ 2 . However, a theoretical determination of s-wave tttt 
scattering lengths which makes use of experimental data has reached a remarkable level of 
precision j§, [Io| : 

m w ai=° = 0.220 ±0.005 ; mW= 2 = -0.0444 ± 0.0010 . (2) 



These values result from the Roy equations [ill, |12J, [13| , which use dispersion theory to relate 



scattering data at high energies to the scattering amplitude near threshold. In a striking 
recent result, this technology has allowed a model-independent determination of the mass 
and width of the resonance with vacuum quantum numbers (the a meson) that appears 
in the tttt scattering amplitude [lij ■ Several low-energy constants of one- loop x~PT are 
critical inputs to the Roy equation analysis. One can take the values of these low-energy 
constants computed with lattice QCD by the MILC collaboration 15, 16[ as inputs to the 



Roy equations, and obtain results for the scattering lengths consistent with the analysis of 
Ref. 0. 

A direct lattice QCD determination of threshold tttt scattering is problematic in two 
respects. First, the occurrence of disconnected diagrams in the 1 = s-wave channel renders 
a determination of that amplitude very costly in terms of computer time, given the current 
state of lattice algorithms, and is thus beyond our current capabilities. As a result, lattice 
QCD efforts have focused on the 1 = 2 channel. The second difficulty is due to the fact 
that lattice QCD calculations are performed on a Euclidean lattice. The Maiani-Testa 
theorem demonstrates that S-matrix elements cannot be determined from lattice calculations 



of n-point Green's functions at infinite volume, except at kinematic thresholds [17|. This 
difficulty was overcome by Luscher, who showed that by computing the energy levels of 
two-particle states in the finite-volume lattice, the 2^2 scattering amplitude can be 



recovered [18|, [19|, |20|, [21|, [22J. The energy levels of the two interacting particles are found 
to deviate from those of two non-interacting particles by an amount that depends on the 
scattering amplitude and varies inversely with the lattice spatial volume. 
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The first lattice calculations of tttt scattering; were performed in quenched QCD [23|, |2 



27|, |2a, |29|, [3d, |3l|, [32j, [33j, |3J, |35|, |36j, |37|, |38j, |39|, |4fl, |4l|, |42( , and the first full-QCD 



calculation of tttt scattering (the scattering length and phase-shift) was carried through by 
the CP-PACS collaboration, who exploited the finite- volume strate gy t o study 1 = 2, s-wave 
scattering with two flavors (n/ = 2) of improved Wilson fermions [43j, with pion masses in 
the range m n ~ 0.5 — 1.1 GeV. The first fully-dynamical calculation of the / = 2 tttt 
scattering length with three flavors (rif = 2 + 1) of light quarks was performed by some 
of the present authors using domain-wall valence quarks on asqt ad- improved staggered sea 
quarks at four pion masses in the range m n ~ 0.3 — 0.5 GeV at a single lattice spacing, 
b ~ 0.125 fm [44|. That work quoted a value of the scattering length extrapolated to the 
physical point of 



-0.0426 ± 0.0006 ± 0.0003 ± 0.0018 



(3) 



where the first uncertainty is statistical, the second is a systematic due to fitting and the 
third uncertainty is due to truncation of the chiral expansion. 

In this paper we update our fully-dynamical mixed-action calculation of the I = 2 tttt 
scattering length. Two recent developments motivate an update: i) we have vastly increased 
statistics at the three light-quark masses studied in the original publication; ii) tttt scatte ring 
has been computed with Mixed- Action X"PT (MAx-PT) at next-to-leading order (NLO) [45j, 
461 ] both for two and three flavors of light quarks. Our updated result is: 



1=2 



-0.04330 ± 0.00042 



(4) 



where the statistical and systematic uncertainties have been combined in quadrature. This 
result is consistent with all previous determinations within uncertainties. 

This paper is organized as follows. In Section [Til details of our mixed-action lattice QCD 
calculation are presented. We refer the reader interested in a more comprehensive treatment 
and discussion to our earlier papers. Discussion of the relevant correlation functions and an 
outline of the methodology and fitting procedures can also be found in this section. The 
results of the lattice calculation and the analysis with two- and three-flavor MA^-PT are 
presented in Section IIHI In Section ITV] the various sources of systematic uncertainty are 
identified and quantified. In Section [V] we conclude. 



II. METHODOLOGY AND DETAILS OF THE LATTICE CALCULATION 



The computation in this paper uses the mixed-action lattice QCD scheme developed by 



LHPC [471. |48J|. Domain- wall fermion propagators were generated from a smeared source on 
rif = 2 + 1 asqtad-improved (49|, [5(| coarse configurations generated with rooted staggered 
sea quarks [Hlj. Hyper cubic-smeared (HYP-smeared) 52, 53, 54, 55] gauge links were used 



in the domain-wall fermion action to improve chiral symmetry (further details about the 
mixed-action scheme can be found in Refs. 56, 5?J )• The mixed-action calculations we 
have performed involved computing the valence-quark propagators using the domain-wall 
formulation of lattice fermions, on each gauge-field configuration of an ensemble of the coarse 
MILC lattices that are generated using the staggered formulation of lattice fermions 58, 5^, 
60 . 61 . 62| and taking the fourth root of the fermion determinant, i.e. domain-wall valence 
quarks on a rooted-staggered sea. In the continuum limit the nj = 2 staggered action has 
an SU (8) l <S> SU(8)r ® U{l)v chiral symmetry due to the four-fold taste degeneracy of each 
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TABLE I: The parameters of the MILC gauge configurations and domain-wall propagators used in 
this work. The subscript I denotes light quark (up and down), and s denotes the strange quark. The 
superscript dwf denotes the bare-quark mass for the domain-wall fermion propagator calculation. 
The last column is the number of configurations times the number of sources per configuration. 

Ensemble b m bm s bmf wf bmff 10 3 x bm res a # of propagators 



2064f21b676m007m050 


0.007 


0.050 


0.0081 


0.081 


1.604 ± 0.038 


468 


X 


16 


2064f21b676m010m050 


0.010 


0.050 


0.0138 


0.081 


1.552 ±0.027 


658 


X 


20 


2064f21b679m020m050 


0.020 


0.050 


0.0313 


0.081 


1.239 ± 0.028 


486 


X 


24 


2064f21b681m030m050 


0.030 


0.050 


0.0478 


0.081 


0.982 ± 0.030 


564 


X 
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"Computed by the LHP collaboration. 



flavor, and each pion has 15 degenerate additional partners. At finite lattice spacing this 
symmetry is broken and the taste multiplets are no longer degenerate, but have splittings 
that are 0(a 2 b 2 ). While there is no proof, there are arguments to suggest that taking 
the fourth root of the fermion determinant recovers the contribution from a single Dirac 
fermion l . The results of this paper assume that the fourth- root trick recovers the correct 
continuum limit of QCD. 

When determining the mass of the valence quarks there is an ambiguity due to the non- 
degeneracy of the 16 staggered bosons associated with each pion. One could choose to 
match to the taste-singlet meson or to any of the mesons that become degenerate in the 
continuum limit. Given that the effective field theory exists to describe such calculations 
at finite lattice spacing, the effects of matching can be described, and removed, by effective 
field theory calculations appropriate to the choice of matching. The quantity b 2 Aj is the 
mass-difference between a valence meson and the staggered taste-singlet meson when the 
valence pion is tuned to be exactly degenerate with the lightest staggered pion. On the 
coarse MILC lattices with b ~ 0.125 fm (and L ~ 2.5 fm) it is numerically determined (in 
lattice units) that 6 2 A 7 = 0.0769(22) [l5|. 

A summary of the lattice parameters and resources used in this work is given in Table [H In 
order to generate large statistics on the existing MILC configurations, multiple propagators 
from sources displaced both temporally and spatially on the lattice were computed. The 
correlators were blocked so that one average correlator per configuration was used in the 
subsequent Jackknife statistical analysis (that will be described later). 

The 7T correlation function, C n (t), and the tttt correlation function C nn {p,t) were com- 
puted, where the number of time slices between the hadronic sink and the hadronic source is 
denoted by t, and p denotes the magnitude of the (equal and opposite) momentum of each 
pion. The single-7r + correlation function is 

<V(*) = £<7r-(*,x)7r + (0,0)) , (5) 

where the summation over x corresponds to summing over all the spatial lattice sites, thereby 
projecting onto the momentum p = state. A 7r + 7r + correlation function that projects onto 



1 For a nice introduction to staggered fermions and the fourth-root trick, see Ref. 63J. For the most 



recent discussions regarding the continuum limit of staggered fermions with the fourth-root trick, see 
Ref. @,S,S,0[6iS0[zi|. 
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the s-wave state in the continuum limit is 



<W(p,t) = £ £e 4 P-( x -^(7r-(t,x)7r-(t,y)7r + (0,O)7r + (0,O)) 

|p|=p x,y 



(6) 



where, in eqs. flSD and (jSJ), vr + (t, x) = tt(t, x)7 5 d(i, x) is an interpolating field (Gaussian- 
smeared) for the 7r + . In the relatively large lattice volumes that we are using, the energy 
difference between the interacting and non-interacting two-meson states is a small fraction 
of the total energy, which is dominated by the masses of the mesons. In order to extract 
this energy difference we formed the ratio of correlation functions, G n + n +(p,t), where 



G m 



C?r+ir+ (g» t) 

C 7T +(t)C 7T +(t) 



CO 

E 

n=0 



-AE„ t 



(7) 



and the arrow denotes the large-time behavior of G n + n + in the absence of boundaries on the 
lattice and becomes an equality in the limit of an infinite number of gauge configurations. In 
G 7T + n +, some of the fluctuations that contribute to both the one- and two-meson correlation 
functions cancel, thereby improving the quality of the extraction of the energy difference 
beyond what we are able to achieve from an analysis of the individual correlation functions. 

The energy eigenvalue E n and its deviation from the sum of the rest masses of the particle, 
AE n , are related to the center-of-mass momentum p n by 



AE n = E n - 2m^ = v l + m\ 



2m„ 



(8) 



In the absence of interactions between the particles, |pcot5| = oo, and the energy levels 
occur at momenta p = 27rj/L, corresponding to single-particle modes in a cubic volume. In 
the interacting theory, once the energy shift has been computed, the real part of the inverse 
scattering amplitude is determined via the Luscher formula [la . 1 191 . l20l . l2lfl . To obtain 
p cot S(p), where S(p) is the phase shift, the magnitude of the center-of-mass momentum, p, 
is extracted from the energy shift, given in eq. (|S]), and inserted into 18, 19, 0, 21, 22]: 



pcot 5(p) 



7T 



L V 2vr J 



(9) 



which is valid below the inelastic threshold. The regulated three-dimensional sum is [22 

1 



|j|<A 

E TT 



J I 2 - V 2 



47rA 



(10) 



where the summation is over all triplets of integers j such that | j | < A and the limit A — > oo 
is implicit. The approximate formula [H, 19, [20, 21] that can be used for L ^> a is 



AEn 



Ana 



m w L 3 



a (a 
l + c lz + C2 ( z 



O 



1 



which relates the ground-state energy shift to the phase shift, with 



Cl 



i IJKA , 

^Et4 



4A 



-2.837297 



TT 



C 2 



^2 Z.^ i:i4 



TT 



(11) 



6.375183 , (12) 
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and a is the scattering length, defined by 



a = lim tang(p) (13) 

For the I = 2 -kit scattering length that we compute here, the difference between the exact 
solution and the approximate solution in eq. ffTTl) is < 1%. For the volumes we consider 
(with L ~ 2.5 fm), the center-of-mass momentum is obviously non-zero and therefore one 
should keep in mind the effective range expansion: 

p cot 5(p) = - + \-rp 2 + OQo 4 ) , (14) 
a — 

where r is the effective range, which appears at O (1/L e ) in eq. (TIT]) , and include the trun- 
cation of eq. CHI) as a source of systematic uncertainty. 



III. DATA ANALYSIS AND CHIRAL AND CONTINUUM EXTRAPOLATION 

A. Results of the Lattice Calculation 

It is convenient to present the results of our lattice calculation in "effective scattering length" 
plots, simple variants of effective-mass plots. The effective energy splitting is formed from 
the ratio of correlation functions 

^w = *(<^Sij) ■ (15) 

which in the limit of an infinite number of gauge configurations would become a constant 
at large times that is equal to the lowest energy of the interacting 7r +, s in the volume. At 
each time-slice, AE n + w +(t) is inserted into eq. (jSJ) (or eq. (fill), to give a scattering length 
at each time slice, a n + n +(t). It is customary to consider the dimensionless quantity given 
by the pion mass times the scattering length, m n a n + n +, where m n (t) is the pion effective 
mass, in order to remove scale-setting uncertainties. For each of the MILC ensembles that 
we analyze, the effective scattering lengths are shown in fig. [TJ The statistical uncertainty 
at each time slice has been generated with the Jackknife procedure. The values of the pion 
masses, decay constants and tttt energy-shifts that we have calculated are shown in Table [Til 



B. Two-Flavor Mixed- Action x-PT at One Loop 

The mixed-action corrections for the I = 2 tttt scattering length have been determined in 
Ref. 45J. It was demonstrated that when the extrapolation formulae for this system are ex- 
pressed in terms of the lattice-physical parameters 2 as computed on the lattice, m n , and /„-, 



We denote quantities that are computed directly from the correlation functions, such clS 777/ 7f . clS lattice- 
physical quantities. These are not extrapolated to the continuum, to infinite-volume or to the physical 
point. 
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FIG. 1: The effective n + n + scattering length times the effective ir mass as a function of time slice 
arising from smeared sinks. The solid black lines and shaded regions are fits with 1-a statistical 
uncertainties tabulated in Table [77J. The dashed lines are estimates of the systematic uncertainty 
due to fitting, also given in Table I77J . 



there are no lattice-spacing-dependent counterterms at 0(b ), 0(b) or 0(m^b ) ~ 0(b ). 
This was explained to be a general feature of the two-meson systems at this order, including 
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TABLE II: The summary table of raw fit quantities required for the two-flavor analysis. The first 
uncertainties are statistical, the second uncertainties are systematic uncertainties due to fitting 
and the third uncertainty, when present, is a comprehensive systematic uncertainty, as discussed 
in the text. 



Quantity 


™ , n 007 

771/ — U.UU ( 


™ , n run 

771/ — U.U1U 


™ , n ooo 

Till — u.uzu 


™ , n o^o 


Fit Range 


8-12 


8-13 


7-13 


9-12 


m v (l.u.) 
h (1-u.) 


0.18454(58)(51) 
0.09273(29) (42) 
1.990(11)(14) 


0.22294(31)(09) 
0.09597(16)(10) 
2.3230(57)(30) 


0.31132(28)(21) 
0.10179(12) (28) 
3.0585(49)(95) 


0.37407(49) (12) 
0.10759(28)(17) 
3.4758(98) (60) 


Fit Range 


11 - 15 


9-15 


10 - 15 


12 - 17 


AE m (l.u.) 

m y= 2 (b ? 0) 

^ 2 (MO) 
5 (b 7^ 0) (degrees) 

\p\/m„ 


0.00779(47) (14) 
-0.1458(78) (25) (14) 
6.1(1.9)(0.7)(0.4) 
-1.71(14)(04) 
0.2032(60)(18) 


0.00745(20)(07) 
-0.2061(49)(17)(20) 
5.23(68)(24)(28) 
-2.181(81)(28) 
0.1836(25)(09) 


0.00678(18)(20) 
-0.3540(68)(89)(35) 
6.53(32) (42) (16) 
-3.01(09)(12) 
0.1480(17) (23) 


0.00627(23)(10) 
-0.465(14) (06) (05) 
6.90(40) (18) (13) 
-3.46(17)(07) 
0.1298(24) (10) 



TABLE III: Summary table for fit quantities extrapolated to the continuum with two-flavor 
MA%PT. The first row corresponds to the overall mixed action correction to the scattering length. 
The uncertainties are discussed in detail in Section HVl The second and third rows are the con- 
tinuum limit scattering length and low-energy constant. The first uncertainties are statistical and 
the second uncertainties are comprehensive systematic uncertainties. 



Quantity 


m t = 0.007 


mi = 0.010 


m t = 0.020 


m t = 0.030 


A (m„al= 2 ) 


0.0033(02)(02)(32)(55) 


0.0030(02)(04)(35)(22) 


0.0023(01)(10)(36)(03) 


0.0018(01)(16)(32)(01) 


mW= 2 (b -> 0) 
fe 2 (6 - 0) 


-0.1491(78)(32) 
5.3(1.9)(1.8) 


-0.2091(49)(34) 
4.83(68) (73) 


-0.356(07)(11) 
6.42(32)(51) 


-0.467(14) (09) 
6.85(40)(27) 



the non-zero momentum states [461 ] . There are additional lattice-spacing corrections due to 
the hairpin interactions present in mixed-action theories, but for our scheme of domain-wall 
valence propagators calculated in the background of the asqtad improved MILC gauge con- 
figurations, these contributions are completely calculable without additional counterterms 
at NLO, as they depend only upon valence meson masses and the staggered taste-identity 
meson mass splitting [45|, |4y] which has been computed [3]. This allows us to precisely 
determine the predicted mixed-action corrections for the scattering lengths at the various 
pion masses used in this work. In two-flavor MA^-PT (i.e. including finite lattice-spacing 
corrections) the chiral expansion of the scattering length at NLO takes the form El 



m n a 



1=2 



(M0) 



m 



1 + 



167T2/2 



3 log 



ml 



A 4 

_ l _ ffU) _ 



where it is understood that and f n are the lattice-physical parameters [46J and 
K) u = - m\ u = 2B ( mj - m u ) + 6 2 A 7 + . . . , 



, (16) 



(17) 



where u denotes a valence quark and j denotes a sea-quark, and we are using isospin- 
symmetric sea and valence quarks, m 



(m uu ) is the mass of a meson composed of two sea 
(valence) quarks of mass rrij (m u ) and the dots denote higher-order corrections to the meson 
masses. Clearly eq. (|T6|) . which contains all (9(m^6 2 ) and 0{b A ) lattice artifacts, reduces to 



S 



the continuum expression for the scattering length in the QCD limit where A 2 ju -> 3 . 
It is worth noting that eq. (fl6]) : and the subsequent expression for the three-flavor theory, 
become the partially-quenched formulae in the continuum limit. Therefore, they are the 
correct extrapolation formulae to use in the case of non-degenerate valence- and sea-quark 
masses, as is implied by eq. ([T5]) and eq. (TTTT) . This modification of the partially-quenched 
formulae can be understood on more general grounds, as mixed-action theories with chirally- 
symmetric valence fermions exhibit many universal features (72| . 



With domain-wall fermion masses tuned to match the staggered Goldstone pion |47l . 148 



one finds A| u = b 2 Aj. The various fit parameters relevant to the two-flavor extrapolation 
are presented in Table [Til For each ensemble we determine a 1 ^ 2 , and then use the chiral 
extrapolation formula to extract a value of the counterterm Z^ 2 (/x = fn), with a statis- 
tical uncertainty determined with the Jackknife procedure. The systematic uncertainties 
are propagated through in quadrature. The results of the two-flavor extrapolation to the 
continuum are shown in Table IHIl 

Fitting to lattice data at the lightest accessible values of the quark masses will optimize 
the convergence of the chiral expansion. While we only have four different quark masses 
in our data set, with pion masses, m n ~ 290 MeV, 350 MeV, 490 MeV and 590 MeV, 
fitting all four data sets and then "pruning" the heaviest data set and refitting provides 
a useful measure of the convergence of the chiral expansion. Hence, in "fit A", we fit the 
/^ 2 (/i = /7r)'s extracted from all four lattice ensembles (m007, mOlO, m020 and m030) to a 
constant, while in "fit B", we fit the l^ 2 (fi = /tiO's from the lightest three lattice ensembles 
(m007, mOlO and m020). In "fit C", we fit the li= 2 (fi = / w )'s from the lightest two lattice 
ensembles (m007 and mOlO). Results are given in Table HVl 



TABLE IV: Results of the fits in two-flavor Mixed- Action %-PT. The values of m n a 1 ^ 2 correspond 
to the extrapolated values at the physical point. The first uncertainty is statistical and the second 
is a comprehensive systematic uncertainty. 



FIT 


&» = U) 


a 1 ^ 2 (extrapolated) 


X 2 /dof 


A 


6.43 ±0.23 ±0.26 


-0.043068 ± 0.000076 ± 0.000085 


1.17 


B 


5.97 ±0.29 ±0.42 


-0.043218 ± 0.00009 ± 0.00014 


0.965 


C 


4.89 ±0.64 ±0.68 


-0.04357 ± 0.00021 ± 0.00022 


0.054 



Taking the range of parameters spanned by fits A-C one finds: 

ll= 2 (fi=U) = 5.4 ±1.4 

m w ai= 2 = -0.04341 ± 0.00046 . (18) 

In Fig. [2] we show the results of our calculation, along with the lowest mass n/ = 2 
point from CP-PACS (not included in our fit). We also show the tree- level prediction and 
the results of our two-flavor fit described in this section. The experimental point shown in 
Fig. [2] is not included in the fit and extrapolation. It is interesting that the lattice data 
indicates little deviation from the tree level ;\PT curve. The significant deviation of the 
extrapolated scattering length from the tree-level result is entirely a consequence of fitting 
to MAxPT at one-loop level. 



3 The counterterm l^ 2 (fi) is, of course, the same counterterm that appears in continuum xPT. 
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m I f 

FIG. 2: m w a 1 ^ 2 vs. m w / f w (ovals) with statistical (dark bars) and systematic (light bars) 
uncertainties. Also shown are the experimental value from Ref. [a]__( diamond) and the lowest 
quark mass result of the nj = 2 dynamical calculation of CP-PACS [43f] (square). The blue band 
corresponds to a weighted fit to the lightest three data points (fit B) using the one-loop MA^-PT 
formula in eq. (|16p (the shaded region corresponds only to the statistical uncertainty). The red 
line is the tree-level x _ PT result. The experimental data is not used in the chiral extrapolation 
fits. 



C. Three-Flavor Mixed- Action x-PT at One Loop 



An important check of the systematic uncertainties involved in the chiral extrapolation is to 
perform the same analysis using three-flavor MA^-PT 45, 4fj| as both the real world and our 
lattice calculation have three active light flavors. In addition to the computations presented 
in Table [IIJ it is necessary to determine masses and decay constants for the kaon and the 
7]. We use the Gell-Mann-Okubo mass-relation among the mesons to determine the rj mass, 
which we do not compute in this lattice calculation due the enormous computer resources 
(beyond what is available to us) required to compute the disconnected contributions. This 
procedure is consistent to the order in the chiral expansion to which we are working. 

The chiral expansion of the 7r + 7r + scattering length in three-flavor mixed-action %PT 



10 



TABLE V: The summary table of quantities required for the three-flavor analysis. A " * " denotes 
that the Gell-Mann-Okubo mass relation among the mesons has been used to determine this 
quantity. The first uncertainties are statistical and the second are systematic (that are discussed 
in the text). 



Quantity 


TO; = 0.007 


m t = 0.010 


mi = 0.020 


mi = 0.030 


Fit Range 


8-14 


9-14 


9-13 


9-13 


m K (l.u.) 
m n (l.u.) * 

nh t /fn * 

rhx/U * 


0.36839(40) (29) 
0.41 182 (43) (36) 
4.447(19)(20) 
5.408(23)(24) 


0.37797(30)(03) 
0.41703(32)(04) 
4.3517(96)(43) 
5.271(11)(05) 


0.40540(31)(32) 
0.43224(33) (46) 
4.246(06) (12) 
5.087(07) (14) 


0.42976(41)(20) 
0.44688(38) (26) 
4.154(11)(05) 
4.927(13)(06) 


S * 

mW-^ 2 (6^0)* 
32(4nfLi= 2 * 


-0.0015(01) 
0.0011(01) 
-0.1470(78) (70) 
6.4(1.9)(1.7) 


-0.0027(00) 
0.0003(01) 
-0.2065(49)(50) 
5.66(67) (68) 


-0.0079(01) 
-0.0012(01) 
-0.353(07) (10) 
7.07(32) (48) 


-0.0130(03) 
-0.0018(01) 
-0.462(14) (08) 
7.44(40) (21) 



takes the form [461: 

7^0) = 



J=2i 



mz 



mz 



**f* 



16tT 2 / 2 



A! 



mz\ . 2 Tl=2 , . 1 



,~2 



31og(-f | -32(4^ 4=^) + -log 



m 



x 



6mt 



- + y 

,4 ^ 



n=l 




9 

(19) 



where m 2 x 



m 2 + 6 2 A/, and 



Mv) 
My) 



■2y 



9(1 -y) 2 
2y 



3(1 -yf 

y 



9(1 -y) A 
y 2 



54(1 -yf 

In addition, it is useful to define the quantities: 

2i\mt 



5(l-y) + (3 + 2y) \n{y) 
l-y)(l + 3y)+y(3 + y) \n(y) 
l-y)(l-7y- I2y 2 ) - 2y 2 (7 + 2y) \n{y) 

(1 - y)(l - 8y - 17y 2 ) - Qy 2 (3 + y) \n(y) 



(47T/ ff ) 




and 



mz 



mz 





(20) 



(21) 



(22) 



87T/2 9 *\f 2 

whose numerical values for the various ensembles are given in Table |Vj 

For the three-flavor analysis, we follow the same procedure of "pruning" the data as in 
the two-flavor analysis, giving the results shown in Table [VTl Taking the range of parameters 
spanned by fits D-F one finds: 

32(4n)L I = 2 ^ = U) 



,1=2 



6.2 ± 1.2 
-0.04330 ± 0.00042 



(23) 
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TABLE VI: Results of the NLO fits in three-flavor Mixed-Action %-PT. The values of m n a[= 2 
correspond to the extrapolated values at the physical point. The first uncertainty is statistical and 
the second is a comprehensive systematic uncertainty. 



FIT 


32(4vr)Li= 2 (^ = / 7r ) 


m v a^ 2 (extrapolated) 


X 2 /dof 


D 


7.09 ± 0.23 ± 0.23 


-0.042992 ± 0.000076 ± 0.000077 


0.969 


E 


6.69 ±0.29 ±0.39 


-0.04312 ± 0.00009 ± 0.00013 


0.803 


F 


5.75 ± 0.63 ± 0.64 


-0.04343 ± 0.00021 ± 0.00021 


0.073 



TABLE VII: Corrections and uncertainties in m^a^ 2 for n/ = 2. 



Quantity 


mi = 0.007 


mi = 0.010 


mi = 0.020 


mi = 0.030 




0.0033(02)(02) 


0.0030(02)(04) 


0.0023(01)(10) 


0.0018(01)(16) 


Afv (™ttg4= 2 ) 


±0.0055 


±0.0022 


±0.0003 


±0.0001 


A mres (m^a 1 ^ 2 ) 


±0.0032 


±0.0035 


±0.0036 


±0.0032 



IV. SYSTEMATIC UNCERTAINTIES 

This section describes the sources of systematic uncertainty that need to be quantified. 



A. Higher-Order Effects in Mixed- Action x-PT 



We rely on the power counting associated with the chiral expansion of the Mixed- Action x?T 
to estimate the size of the lattice-spacing artifacts arising at 0{m\b 2 ). To be conservative, 
we have estimated these corrections to be of the general size 



0{m%b 



2-irmi 



b 2 A 1 



(4vr/ 7r ) 4 (47TA) 



(24) 



We treat these estimates as uncertainties in the predicted NLO MA^PT corrections which 
can be determined from eq. ( JT6"1) and eq. (TT9l . We provide these predicted corrections and 
their uncertainties in the form 



A 



MA 



1=2 



MA 



1=2 



(25) 



The values of these corrections are shown in Tables IVHI and IVIIII The first uncertainty 
in these corrections is statistical and is associated with the meson masses, decay constants 
and the taste-identity mass splitting, b 2 Ai. The second uncertainty is the power counting 
estimate of the higher-order corrections of 0(m^.b 2 ) as estimated in eq. (j2"4l . The calcula- 
ble corrections to m^a 1 ^ 2 at C(m 2 6 2 ,6 4 ) are 2.3%, 1.5%, 0.65% and 0.39% effects for the 
007, 010, 020 and 030 ensembles, respectively, from which we conclude that the C(m 4 6 2 ) 
contributions are significantly less than ~ 1%. 
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TABLE VIII: Corrections and uncertainties in m^a^ 2 for rif = 2 + 1. 



Quantity 


mi = 0.007 


mi = 0.010 


mi = 0.020 


mi = 0.030 




0.0012(01)(02) 


0.0004(01)(04) 


-0.0015(03)(10) 


-0.0027(05)(16) 


A FV (m w 4= 2 ) 


±0.0024 


±0.0005 


±0.0001 


±0.00006 


A mres (m,a^ 2 ) 


±0.0032 


±0.0035 


±0.0036 


±0.0032 



B. Finite- Volume Effects in Mixed- Action x-PT 



The universal relation between the two-particle energy levels in a finite volume and their 
infinite-volume scattering parameters receives non-universal corrections which are exponen- 
tially suppressed by the lattice size and dominated by the lightest particle in the spectrum. 
These scale generically as e~ m ^ L 7a, \74 \ . In Ref. [75(, the leading exponential volume cor- 
rections to p cot 5(p) were determined in the I = 2 tttt scattering channel in %PT. However, 
in order to determine the leading finite-volume corrections to this mixed-action calculation, 
hairpin diagrams present in the mixed-action theory must also be included. For the 1 = 2 tttt 
system, there are additional hairpin diagrams present in the t and u channel scattering di- 
agrams jisl ]. The finite- volume corrections from these diagrams are larger than those in 
continuum x?T, but are opposite in sign and therefore the overall magnitude of the correc- 
tion is similar to that given in Ref. [75j. We note that as these contributions vanish in the 
continuum limit, they are actually finite- volume finite- lattice-spacing corrections, and not 
just finite-volume corrections, and hence scale as 6 2 exp(— m n L) at small lattice spacing. 

As with the mixed-action lattice-spacing corrections, we denote these finite-volume mod- 
ifications as 



A 



FV 



1=2 



FV 



1=2 



ooV 



(26) 



and they are shown in Tables IVlTl and IVIII1 However, one should take note that the effective- 
range contribution to pcot 5(p), which behaves as a power-law in the lattice size (and there- 
fore is parametrically enhanced over the exponential corrections) is not included in the 
extraction of the scattering lengths. While the exponential modifications are numerically 
larger than our estimate of the effective-range contributions at the light pion masses (see 
below), the values of Api/(m T a^ 2 ) shown in Tables IVHl and IVIIII are used as estimates of 
the uncertainties due to higher-order finite-volume effects. 



C. Residual Chiral Symmetry Breaking 

The mixed-action formulae describing tttt scattering determined in Refs. EE ra have as- 
sumed that the valence fermions have exact chiral symmetry, up to the quark-mass correc- 
tions. The domain-wall propagators used in this work have a finite fifth-dimensional extent 
and therefore residual chiral symmetry breaking arising from the overlap of the left- and 
right-handed quark fields bound to the opposite domain walls. Due to the nature of this 
residual chiral symmetry breaking in the domain-wall action, the leading contributions can 



be parameterized as an additive shift to the valence-quark masses [6d . |62j , 

mf w f — > mf w f + m res . (27) 
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A full treatment of these effects involves three new spurion fields in the effective field the- 
ory [76] but this is not necessary for estimating the size of these contributions to the nir 
scattering lengths. By expressing the calculated scattering lengths and extrapolation for- 
mulae in terms of the lattice-physical meson masses and decay constants, the dominant 
contributions from residual chiral symmetry breaking are included, leaving corrections at 
higher orders in the chiral expansion. There will be new operators similar to the Gasser- 
Leutwyler operators 77J in the chiral Lagrangian, for example 



C = 2B U str (0 M £0"£t) str (m res £ t + £mj es ) 



L 6 str (m^ + EmJ 



str [m res I^ + £mj es ) + 



(28) 



Naive dimensional analysis [78j can be used to estimate the size of the corrections due to 
these new operators, which in the case of the I = 2 tttt system are given by 

,4 



I=2\ 



inmi m r 



(29) 



There will be additional operators with two insertions of m res in the place of m q , but these 
are < 20% of the uncertainty already estimated for the residual chiral symmetry breaking. 
These uncertainties are denoted by 



A mrea (m n ai n 2 ) = m T d 



1=2 



.1=2 



(30) 



and are shown in Tables IVHI and IVIIII 



D. Two Loops Effects 



The two-loop expression for the scattering length 0, @] is given, in the continuum limit of 
QCD, by 



.1=2 



ml 



ml 



8^ 



167T 2 /. 2 



3 log 



ml 



ml 



64tt 4 /, 4 



,2\ 2 



i - C>) 



3 4 log 5) + ev)iog^ + e^) 



6 



I* 



(31) 



where it® and l^J are linear combinations of undetermined constants that appear in the 
0(p 4 ) and 0(p 6 ) chiral Lagrangians [2j,|4|. Fitting all four data points allows for an extraction 
of the three counterterms with x 2 /dof = 0.26. From the 68% confidence-interval error 
ellipsoid we find an extrapolated value of: 

(32) 



m-r: a 



1=2 



-0.0442 ± 0.0030 . 



While it is gratifying to have a determination of the scattering length at two-loop level that 
is consistent with the one-loop result, there are several caveats: i) the two-loop expression 
in MA^-PT does not yet exist and therefore the determination in eq. ( |3T1) contains lattice- 
spacing artifacts at lower orders in the chiral expansion than in the one- loop result; ii) This 
value is clearly strongly dependent on the heaviest quark mass, which is, at best, at the 
boundary of the range of validity of the chiral expansion. A reliable two-loop determination 
will have to await further lattice data at quark masses closer to the chiral limit than we 
currently possess. 
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E. Range Corrections 

It is straightforward to show that the range corrections enter at (9 (L -6 ) in eq. (ITT1) . As- 
suming that the effective range is of order the scattering length (the scattering length is of 
natural size), we expect a fractional uncertainty of {m^a) 2 p 2 /2m\ due to the omission of 
range corrections. For the ensembles that we consider, this translates into an 0.5% uncer- 
tainty in m n a I n ^ 2 . Allowing for the effective range to exceed its natural value by a factor of 
two, we assign a 1% systematic uncertainty to m^a 1 ^ 2 determined on each ensemble. 

F. Isospin Violation 

The calculation we have performed is in the limit of exact isospin symmetry, as are the 
extrapolation formula we have used to analyze the results. The conventional discussion of 
the scattering length is in the unphysical theory with e = and m u = = m, with 
77^ = = 139.57018 ± 0.00035 MeV and f n = f n + = 130.7 ± 0.14 ± 0.37 MeV. Hence 
m w + 1 f n + = 1.0679 ± 0.0032, where the statistical and systematic uncertainties have been 
combined in quadrature. We extrapolate the results of our lattice calculations to this value. 

Unfortunately, we are presently unable to make precise predictions for the real world in 
which isospin breaking occurs at the few-percent level. Extrapolation to the isospin-averaged 
pion mass (as opposed to the charged pion mass), would introduce a shift of ~ 2% in m^a 1 ^ 2 . 
This is larger than the uncertainty we have determined at the charged pion mass. It is clear 
that in order to make predictions for real-world quantities at the ~ 1% level from lattice 
QCD calculations, isosp in-breaking and electromagnetism will need to be incorporated into 
the lattice calculation. 

V. DISCUSSION 

We have presented results of a lattice QCD calculation of the I = 2 scattering length 
performed with domain-wall valence quarks on asqt ad- improved MILC configurations with 
2+1 dynamical staggered quarks. The calculations were performed at a single lattice spacing 
of b ~ 0.125 fm and at a single lattice spatial size of L ~ 2.5 fm with four values of the light 
quark masses, corresponding to pion masses of ~ 290, 350, 490 MeV and 590 MeV. High 
statistics were generated by computing up to twenty-four propagators per MILC configura- 
tion at spatially- and temporally-displaced sources. We used one-loop MA^-PT with two and 
three flavors of light quarks to perform the chiral and continuum extrapolations. Our predic- 
tion for the physical value of the I = 2 nn scattering length is m^a 1 ^ 2 = — 0.04330±0. 00042, 
which agrees within uncertainties with the (non-lattice) determination of CGL but we 
emphasize once again that our result rests on the assumption that the fourth-root trick 
recovers the correct continuum limit of QCD. In Table IIXI and fig. [3] we offer a comparison 
of our prediction with other determinations. What has enabled such an improvement in 
precision over our previous result on the coarse MILC lattices is the recent understanding 
of the lattice-spacing artifacts accomplished with mixed-action chiral perturbation theory. 

While it will be quite useful to have results at another lattice spacing and at another 
lattice volume, we have reached the level of precision where we require knowledge of isospin 
violating effects in order to further reduce the uncertainty in the physical tttt scattering 
lengths; i.e. those that can be compared to experiment. One somewhat surprising result of 
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TABLE IX: A compilation of the various calculations and predictions for the 1 = 2 nn scattering 
length. The prediction made in this paper is labeled NPLQCD (2007). Also included are the 
experimental value from Ref. f^J (E 865 (2003)), the previous determination by NPLQCD fJJ] 
(NPLQCD (2005)), two indirect lattice results from MILC 0, 0/ (the stars on the MILC results 



indicate that these are not lattice calculations of the 1 = 2 tttt scattering length but rather a hybrid 
prediction which uses MILC's determination of various low-energy constants together with the Roy 
equations), and the Roy equation determination of Ref. /£/ (CGL (2001)). 



m T ai n 2 



XPT (Tree Level) -0.04438 
NPLQCD (2007) -0.04330 ± 0.00042 

E 865 (2003) -0.0454 ± 0.0031 ± 0.0010 ± 0.0008 

NPLQCD (2005) -0.0426 ± 0.0006 ± 0.0003 ± 0.0018 
MILC (2006)* -0.0432 ± 0.0006 

MILC (2004)* -0.0433 ± 0.0009 

CGL (2001) -0.0444 ± 0.0010 



our analysis is that one of the dominant sources of systematic uncertainty in our calculation 
is due to residual chiral symmetry breaking in the domain-wall valence quarks for the lattice 
parameters we have chosen. Clearly this systematic can be reduced by improving our choice 
of domain-wall parameters. 

Lattice QCD is currently in a precision age insofar as single-particle properties are con- 
cerned. The precise prediction for the intrinsic two-particle property presented here is a 
remarkable demonstration of the power of combining a lattice QCD calculation with the 
model-independent constraints of chiral perturbation theory. 
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FIG. 3: Bar chart of the various determinations of the 1 = 2 rnr scattering length tabulated 
in Table IIXI We reiterate that the stars on the MILC results indicate that these are not lattice 
calculations of the 1 = 2 tttt scattering length but rather a hybrid prediction which uses MILC's 
determination of various low-energy constants together with the Roy equations. 
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